Analytical solutions for the Rabi model 
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The Rabi model that describes the fundamental interaction between a two-level system with a 
quantized harmonic oscillator is one of the simplest and most ubiquitous models in modern physics. 
However, this model has not been solved exactly because it is hard to find a second conserved 
quantity besides the energy. Here we present a unitary transformation to map this unsolvable Rabi 
model into a solvable Jaynes-Cummings-like model by choosing a proper variation parameter. As a 
result, the analytical energy spectrums and wavefunctions including both the ground and the excited 
states can be obtained easily. Moreover, these explicit results agree well with the direct numerical 
simulations in a wide range of the experimental parameters. In addition, based on our obtained 
energy spectrums, the recent experimental observation of Bloch-Siegert in the circuit quantum 
electrodynamics with the ultrastrong coupling can be explained perfectly. Our results have the 
potential application in the solid-state quantum information processing. 

PACS numbers: 42.50.Pq 



I. INTRODUCTION 

The Rabi model, which was introduced over 70 years 
ago, describes the important interaction between a two- 
level system and a quantized harmonic oscillator (bosonic 
mode) [l] . The corresponding Hamiltonian of this model 
is written as 



Hr = cua [ a + — ila z + ga x (a' 



(1) 



where a^ and a are creation and annihilation operators 
for the quantized harmonic oscillator with frequency u>, 
Oi (i = x, y, z) are the Pauli spin operators, VL is the res- 
onant frequency between two levels, and g is the interac- 
tion strength. In modern physics ranging from quantum 
optics H, condensed-matter physics [3j to quantum infor- 
mation [J| , the Rabi model has been one of the simplest 
and most ubiquitous models. Despite its old age and 
central importance, the Rabi model has not been solved 
exactly because it is hard to find a second conserved 
quantity besides the energy. The well-known method 
to obtain the energy spectrums and the wavefunctions 
for the Rabi model is the numerical diagonalization in a 
truncated finite-dimensional Hilbert space [5|. However, 
these numerical results are difficult to extract the funda- 
mental physics of the Rabi model |6H12| and to precisely 
control the experimental parameters to process quantum 
information |13| . 

To overcome the problem in the analytical considera- 
tions, the rotating-wave approximation, which is valid in 
the regime g <C u, was proposed to rewrite Hamiltonian 



© as [3 

H 3C 



uja^a - 



-£la z + g(cr_(v + <r + a). 



(2) 
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Interestingly, the Jaynes-Cummings model with U(l) 
symmetry has a conserved quantity C = a^a + \{u z + 
1), and thus, it can be solved easily in the subspaces 
{\\, n) , \l,n + 1)}. Meanwhile, the Jaynes-Cummings 
model can successfully describe quantum dynamics of op- 
tical cavity electrodynamics with strong coupling. How- 
ever, in the recent investigations about the solid-state 
quantum electrodynamics with the ultrastrong coupling 
(g ~ O.lw) [15B24J. the rotating-wave approximation 
breaks down and the system's dynamics must be gov- 
erned by the Rabi model. Irish put forward a generalized 
rotating-wave approximation to solve the Rabi model. 
This method can be used to derive the analytical energy 
spectrums of Hamiltonian ([1} in the ultrastrong coupling 
successfully [25|]. However, this method works reason- 
able only for the negative detuning (Q < w) |2a - [28| . 
Recently, Braak used the property of Z2 symmetry of 
the Rabi model to obtain its analytical solutions, which 
are, however, dependent of the composite transcendental 
function defining through its power series in the interact- 
ing strength g and the frequency u> [29( . In our previous 
work, we put forward a generalized variational method to 
analytically obtain the ground-state energy of the Rabi 
model, which agrees well with the numerical simulation 
in all regions of the resonant frequency of the two-level 
system including the negative detuning (Q < w), the res- 
onant case (Q = w), and especially the positive detuning 
(Q > lu). Unfortunately, our introduced method cannot 
be used to consider the excited-state energy spectrum 
30]. 

In this paper, we present a simple and straightforward 
method to solve the Rabi model in both ground and ex- 
cited states. We first map the unsolvable Rabi model ([T]) 



into a solvable Jaynes-Cummings-like model by choos- 
ing a proper variation parameter, as shown in Section 
II. Thus, the analytical energy spectrums including the 
ground and the excited states are obtained in section III. 
These derived analytical results agree well with the direct 
numerical simulations in a wide range of the experimental 
parameters. Moreover, the recent experimental observa- 
tion of Bloch-Siegert [2C)||, which is just the energy shift of 
the level transition, in the circuit quantum electrodynam- 
ics with the ultrastrong coupling can be well explained. 
In section IV, the analytical wavefunctions and the cor- 
responding experimentally-measurable physics quantities 
such as the mean photon number are derived. Also, the 
obtained mean photon number can agree well with the 
numerical simulation. Finally, some conclusions are re- 
marked in Section V. 



terms for a and a^ in Eqs. ([9]) and (fTO)) can be eliminated. 
As a result, the effective Hamiltonian ([5]) reduces to the 
form 

He = uia'a — (Xuj + g)a z (a* + a) + uiX +2Xg + 

^Q{a x G (N) + ia y [Fi(N)a) - aF x (N)]} .(11) 

In the eigenstates of a x with a x |±a;) = ± |±x), the 
Pauli spin operators become u z — > — (r + + r_) and 
o y — > — i(r+ - r_), where t z = \+x) (+x\ — \—x) (—x\, 
t + = \+x) (—x\ and t_ = |— x) (+x|, and the effective 
Hamiltonian (jll[) turns into 

H E = uja^a + uuX 2 + 2Xg + -flG (n)T z + 

R r (r + a + T-a)) + R ar {T + a) + r-a), (12) 



II. MAPPING INTO THE 
JAYNES-CUMMINGS-LIKE MODEL 

When performing a rotation around the y axis, the 
Rabi model becomes 



where 



H R = ujc^a + —a x - g(a) + a)a z 

Under a unitary transformation 

U = cxp[Aer 2 (a' — a)] 



(3) 



(4) 



with A being the dimensionless parameter determined 
by the following calculations, an effective Hamiltonian 
is given by He = UHW , namely, 



where 



He — Hi + H-2. + H3 , 



H\ = ua'a — Xijja z (a* + a) + ojX 



H 2 = -g[<r z (<J +a)-2X], 



(5) 

(6) 

(7) 



O 



H 3 = — {a 2; cosh[2A(a t -a)]+icr a sinh[2A(a t -a)]}. (8) 

Since cosh(y) and sinh(y) are the even and odd 
functions respectively, the terms cosh [2 A (a' — a)] and 
sinh[2A(at — a)] can be expanded as 

cosh[2A(a t -a)] = G (A r ) + Gi(N)a^ + a 2 Gi{N) + . .(9,) 
sinh[2A(a t -a)]=Fi(A)a t -aFi(iV) + 

F 2 (A0a t2 -a 2 F 2 (A0+..., (10) 

where G,(A) (i = 0,1,2,- ••) and Fj(N)(j = 1,2,- ■•) 
with N — a) a are the coefficients that depend on the 
dimensionless parameter A and the photon number n. In 
general, the multi-photon process is weak in the Rabi 
model [2]. It means that the terms of the high-order 



G (n) = (n|cosh[2A(a t -a)]|n), (13) 



R r = Xoj + g - -fi/i(m, n), 



R ar = Xuj + g+ -flfi (m, n) 



(14) 



(15) 



m\Fi(N)ai \n)/y/n + l 



with fi(m,n) = 

(n + l\ Fi(N)(v \n) /-y/n+T. It is straightforward to cal- 
culate that G Q (n) = /n|e 2A ( Qt - Q ) +e- 2A ( Qt - Q )|n\/2 = 



n\e 



2A(a t -a)l 



"L„(4A 2 ) and /i(m,n) 



fi(n + l,n) = (n + l|sinh[2A(at-a)|n) 



2A 



.-2A 2 rl 



L„(4A )/(n + 1), where L n (y) is the La- 



guerre polynomial and L\(y) is the associated Laguerre 
polynomial. Thus, if the dimensionless parameter A is 
chosen as Xuj + g + if2/i(n + 1, n) = 0, namely, 



(Xuj + g) + 



nx 

n + l 



-2\ z t\ 



L' n (^ 2 



0. 



(16) 



the effective Hamiltonian (|12[) becomes a Jaynes- 
Cummings-like model 

H E = uja)a + ujX 2 + 2Xg+ -QG Q (n)T z + R r (T + a + r^a^) 7 

(17) 
where R r = 2(Xuj + g). By means of Hamiltonian (|17[) . 
the analytical energy spectrums and wavefunctions for 
the Rabi model can be obtained easily in the subspaces 
{|-hr,n) , \—x,n + 1)} . 

Before proceeding, the solution of the dimensionless 
parameter A in Eq. (fir?)) should be analyzed since it has 
a crucial role in obtaining the explicit energy spectrums 
and wavefunctions. In general, the nonlinear equation 
(1161) cannot be solved analytically. However, in the cur- 
rent experimental setup with the ultrastrong coupling 




parameter becomes A ~ —g/{uj + f_). 
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FIG. 1: (Color online) The energy spectrums of the ground 
and excited states for the Rabi model as a function of the 
interaction strength g for Q = 0.5w (a), Q = I.Olj (b), and 
_. — 1.5uu (c). The lowest line reflects the ground-state energy 
spectrum. In all figures, the black solid lines represent the ex- 
act numerical results, while the red open symbol corresponds 
to the analytical results obtained in the paper. 



(g < 0.5a;), the numerical result shows that the dimen- 
sionless parameter A is small compared with the unit. 
Thus, the associated Laguerre polynomial is given by 
L^(4A 2 ) ~ n + 1 since L\(y) = n + 1 + J2j>o c jV J i an d 

Eq. ([IB]) becomes (Xuj + g) + f_Ae~ 2A = 0, which leads 
to a solution 



A 



nexp[-2(^)>] 



(18) 
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FIG. 2: (Color online) The energy spectrums of the ground 
and excited states for the Rabi model as a function of the 
resonant frequency f. for g = O.luu (a) and g — 0.3a; (b). The 
lowest line reflects the ground-state energy spectrum. In all 
figures, the black solid lines represent the exact numerical re- 
sults, while the red open symbol corresponds to the analytical 
results obtained in the paper. 

III. ANALYTICAL ENERGY SPECTRUMS 

In terms of the Jaynes-Cummings-likc Hamiltonian 
(|17[l and Eq. (fl8|) . the ground-state energy spectrum 
can be written as 



E G = luX 2 + 2Xg - -fie 



-2\ z 



(19) 



For the weak interaction strength g, the dimensionless and the excited-state energy spectrum can be given by 
I 



E±, n = (n + ^)oj + oj\ 2 + 2\g+^e- 2x2 [L n (4\ 2 )-L n+1 (4\ 2 )} 



±-J{ U ^—[L n {W) + Wi(4A 2 )]} 2 + 4[(Aw + g)V^+T -^= 

2y 2 V" + 1 



LU^w 



r 



(20) 



The ground-state energy spectrum in Eq. (|19p is iden- 
tical to the result derived from the generalized variational 
method [30(, which is, however, invalid for calculating the 
excited-state energy spectrum. In addition, for the weak 
resonant frequency f_, Eqs. (flT)]) and (j_0"]) reduce to the 
results obtained by means of the generalized rotating- 
wave approximation 25]. In Fig. 1, the energy spectrums 
of the ground and excited states for the Rabi model as a 
function of the interaction strength g for f. = 0.5a; (a), 
f. = l.Ow (b) and SI = 1.5a; (c) are plotted, compared 
the explicit results in Eqs. (flT))) and (|2U|) with the di- 
rect numerical simulation. This figure shows that our 
analytical energy spectrums including both the ground- 



and excited states agree perfectly with the direct numer- 
ical simulation in the current experimental setup with 
the ultrastrong coupling (g < 0.5a;). Moreover, these re- 
sults arc valid for all parameter regimes with the negative 
detuning (f. < ui), the resonant case (i. = ui), and espe- 
cially the positive detuning (f_ > w). This conclusion can 
be also drawn from Fig. 2, in which the energy spectrums 
as a function of the resonant frequency f2 for g = 0.1a; 
(a) and g — 0.3a; (b) are plotted. For the generalized 
rotating-wave approximation, the derived energy spec- 
trums can agree well with the numerical simulation in the 
case of the negative detuning (J. < u>). However, with 
the increasing of the resonant frequency CI, this method 
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FIG. 3: (Color online) The Bloch-Siegert (BS) shift as a func- 
tion of the interaction strength g with uj/2-k — 8.13 GHz and 
0/27T = 4.25 GHz. In this figure, the black solid lines repre- 
sent the exact numerical results, while the red open symbol 
corresponds to the analytical results obtained in the paper. 



breaks down and the error is increased linearly. 

In addition, by means of Eqs. (fT9| and ([20)) . the Bloch- 
Siegert shift, which is just the energy shift of the level 
transition, can be also obtained. In Fig. (3), the small- 
est Bloch-Siegert shift with the transition E-,o — > Eg as 
a function of the interaction strength g is plotted. The 
other parameters are the same as those in the recent ex- 
periment of circuit quantum electrodynamics with the 
ultrastrong coupling [20], namely, lj/2tt = 8.13 GHz and 
0/27T = 4.25 GHz. If g/w = 0.1, the smallest Bloch- 
Siegert shift is 50 MHzjwhich agrees well with the ex- 
perimental observation [20]. 



IV. ANALYTICAL WAVEFUNCTIONS AND 
MEAN PHOTON NUMBER 

The wavefunctions of the Rabi model {l} for the 
ground and excited states also need to be discussed. For 
the ground state, the corresponding wavefunction can be 
evaluated immediately as 



\fo, 



= F 4 "v, 



-\<7 Z (a' —a) 



\0,-x) 



(21) 



For the excited state, the wavefunction for the Jaynes- 
Cummings-like Hamiltonian (| 17[) is given by 



\+,n) = cos6 n \n, +x) + sm9 n \n + 1, —x) 
\ — ,n) = — sin^n \n, +x) + cos6 n \n + 1, —x) 



where 



1 _ : 2R, 

2 tan_ E— 



E_ 



+,n 



(22) 



(23) 
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FIG. 4: (Color online) The mean photon numbers of the 
ground and excited states for the Rabi model as a function 
of the resonant frequency SI for g — O.Icj. The lowest line 
reflects the ground-state mean photon number. In this figure, 
the black solid lines represent the exact numerical results, 
while the red open symbol corresponds to the analytical re- 
sults obtained in the paper. 



where E±^ n are given by Eq. (I2TJ)) . Thus, the excited- 
state wavefunction for the Rabi model is given by 



\Ve) 



e ^ ye -\a z (a^a) |_ ) m _ 1) } ( m = 1 ) 3 j . . . ) 

e T L ^y e - A ' T -( Qt - a ) |+, m - 1) , (m = 2, 4, • • • ) ' 

(24) 
where m is the number of the excited state. 

Based on the wavefunctions in Eqs. (|21|) and (|24j). 
the experimentally measurable physics quantities can 
be well derived. For example, the ground-state mean 
photon number can be given by (a^a) = (</?o| a^a \fo) = 

(0, -x\ e A M at - a ) e T IL<T "afaeT CT « e - A ^( at - Q ) |0, -x), 
namely, 



<« f «>o = ^ 



(25) 



For the excited state, the mean photon number is given 
by 



<at a ) £ 



A* + A 2 + 



l+2Atan9„ 



v + a 2 + i^^A 



l+tan 2 6„ 



-, (m = 1,3,---) 

^E (m = 2,4,-..) 
(26) 
where /i = (m — l)/2 and v = (m — 2)/2. 

In Fig. 3, the mean photon numbers of the ground and 
excited states for the Rabi model as a function of the res- 
onant frequency $1 for g = O.lw are plotted. This figure 
shows that the analytical results in Eqs. (|2"5)l and (|2"6")l 
agree well with the numerical simulation. It implies again 
that the effective Hamiltonian (|T7l) can describe the cur- 
rent experimental setup with the ultrastrong coupling. 
Eq. (|2"5j) also shows that the ground-state mean photon 
number for the Rabi model depends on all parameters 



including the frequency u> of the quantized harmonic os- 
cillator, the interaction strength <?, and especially, the 
resonant frequency Q. It is quite different from the re- 
sult derived from the generalized rotating-wave approx- 
imation that the ground-state mean photon number is 
independent of the resonant frequency f2 [25J . 



V. CONCLUSIONS 

In summary, we have presented a unitary transforma- 
tion to map the unsolvable Rabi model into a solvable 
Jaynes-Cummings-like model in the dress-state represen- 
tation. As a result, the analytical energy spectrums and 
wavefunctions including both the ground and the excited 
states can be obtained easily. Moreover, our results agree 
perfectly with the direct numerical simulations in a wide 



range of the experimental parameters and are valid for 
all regions of the resonant frequency of the two-level sys- 
tem including the negative detuning (f2 < ui), the reso- 
nant case (Q = uj), and especially the positive detuning 
(51 > lo). Our results can also explain the recent experi- 
mental observation in the circuit quantum electrodynam- 
ics with the ultrastrong coupling. 
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